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^ . Within the Standard Model, we investigate the semi-leptonic weak decays of J/V'- The 

' various form factors of J/tp transiting to a single charmed meson (-D|^'^j) are studied in the 

00 ' framework of the QCD sum rules. These form factors fully determine the rates of the weak 

semi-leptonic decays of J /if: and provide valuable information about the non-perturbative 
(~| . QCD effects. Our results indicate that the decay rate of the semi-leptonic weak decay mode 

J/^ Di*^^ +e+ + Ue is at order of IQ-^". 

^: 
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r-| ; PACS numbers: 13.20.Gd, 13.25.Gv, 11.55.Hx 
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> ' I. INTRODUCTION 

oo ■ 
cn 

Although strong and electromagnetic decays of J ftp have been extensively studied for several decades, 
I both experimental and theoretical investigations of weak decays of J/ip are much behind. Due to smallness 
. of the strength of weak interaction, the weak decays of the J/tp are rare processes. Sanchis-Lonzano 
suggested to search for these rare decays whose sum of branching ratios were estimated to be at the order 



^ ■ of 10 * Such processes hardly drew much attention because the database was far from reaching such 

H ; 

5^ , accuracy. Thus, for a long time, few further researches on this topic were done. Thanks to the progress of 
accelerator and detector techniques, more accurate measurements may be carried out, thus the interest 
on weak decays of J/tp has been revived. The BES collaboration indeed starts to measure some rare 
weak decays of J/^p and eventually sets an upper bound on the branching ratio of J/^p D + e + Ve at 
order of 10~^ by using 5.8 X 10"^ J/ij datab ase [4]. The forthcoming upgraded BESIII can accumulate 
10^*^ J/ip per year Q], which makes it marginally possible to measure such weak decays of J/ip, at least 
one may expect to observe not-null such events. Thus, more careful theoretical investigation on these 
decays seems necessary. 

Indeed, the weak decays of heavy quarkonium like J/ip offer an ideal opportunity of studying non- 
perturbative QCD effects, because such systems contain two heavy constituents of the same flavor. The 
situation is quite different from that for heavy mesons which contain only one heavy constituent, and 
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the non-perturbative effects miglit be attributed to the hght fiavor, thus the heavy quark effective theory 
(HQET) applies. Moreover, for the weak decay of a vector meson, the polarization effect may play a role 
to probe the underlying dynamics and hadron structure 

The weak decay of J/ip is realized via the spectator mechanism that the charm quark (antiquark) 
decays and the antiquark (quark) acts as a spectator. The characteristic of the decay modes is that 
the final state contains a single charmed hadron. The theory of weak interactions has been thoroughly 
investigated and the effective hamiltonian at the quark level is perfectly formulated. The main job of 
calculating the rates of the semi-leptonic decays of J is to properly evaluate the hadronic matrix 
elements for J/'ijj D^*\ namely the transition form factors which are obviously governed by non- 
perturbative QCD effects. The main aim of this work is to calculate the J/V' — > ^\ds) form factors in 
the QCD sum rules. 

The weak decay of heavy quarkonium has been studied by virtue of heavy quark spin symmetry 
In that framework, the transition form factors of a heavy quarkonium to heavy pseudoscalar and 
vector mesons are parameterized by a universal function 7712 (^'i •V2) in analog to the Isgur-Wise function 
for the heavy meson transitions. However, the non-recoil approximation 7/12(^1 • V2) ~ 1 was used in 
Ref.[l|, which would bring up uncontrollable uncertainties to the estimation of decay widths. It seems 
helpful to re-investigate these processes based on a more rigorous theoretical framework. Motivated by 
the arguments, in this work we will calculate the form factors for heavy quarkonium J /ip decays into a 
pseudoscalar or vector meson in the QCD sum rules. 

As a matter of fact, many authors have tried to evaluate the transition form factors for the heavy meson 
and quarkonium system in various approaches, such as the sirnple quark model {4], light-front approach 

the QCD sum rules Ql) the perturbative QCD approach [s] and etc. The QCD sum-rule approach, 
which is rooted in the quantum field theory and fully relativistic, is considered to be one of the effective 
tools for analyzing hadronic processes Q]. Besides evaluation of hadron spectra, the QCD sum-rule 
technique has been applied to calculate the pion electromagnetic form factor at intermediate momentum 
transfer [s, [lO|, various weak decay channels [h, [l^, the coupling constant of the strong interaction 
[1^ and even to determine the light cone distribution amplitudes of hadrons [l^ . The advantage of this 
method is that the non-perturbative QCD effects are included in a few parameters such as the quark- 
and gluon-condensates which have evident physical meaning jisl . 

After this introduction, we will firstly display the effective Hamiltonian relevant to the semi-leptonic 
decays of J /iJj to ) and the sum rules for form factors in section [Hi The Wilson coefficients of 

various operators which manifest the perturbative QCD effects are also calculated in this section with 
the help of operator product expansion (OPE) technique. The numerical analysis on the form factors are 
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performed in section [Till The decay rates of semi-leptonic decay J/V' -^d(s) ^ comparison of 

our results with that obtained based on other approaches are presented in section IIVI In the last section 
we draw our conclusion. 



II. j/iA ^ d'-*^1^ transition form factors in the qcd sum rules 

A. Definitions of J/tJj ^d{s) transition form factors 



For the semi-leptonic decays J/ip ^ D^^^^^l~^vi, the effective weak Hamiltonian is given by 



-U,fj{c ^ s{d)l9i) = - 75)cz7;7^(l - 75)/. 



(1) 



In order to calculate the rate of a semi-leptonic decay, the essential ingredient is the hadronic matrix 



element (D^^^^|s7^(l — 75)c| J/?/') which is parameterized by various form factors 



(-Dd(s)(p2)k7A«(l - 75)c| J/V'(e,Pi)) 

2V{q^) 



+ i{m^ + uid) 



e • q 



+ i- 



2 ^2 
2 ''?/^ 



(^d(^)(e2,P2)k7M(l - 75)c| J/V'(ei,pi)) 



+ ^2 _^^2 W/?KP2[-^3(g')eie2 • q - ^4(g')erei • <l] 



+ 



ml - ml, ■ 

+ 2 — qf^^^iq ) 



(2) 



(3) 



where the convention Tr[7^7^7p7o-75] = ^ie^upa is adopted. For a transition of J/ip into a charmed 
pseudoscalar meson which is induced by the weak current, there are four independent form factors: 
V, Aq, Ai, A2; while there are ten form factors for J/?/^ transiting to a charmed vector meson which 
are parameterized as Ai{i = 1,2,3,4), Vj{j = 1,2,3,4,5,6). It is worthwhile to emphasize that the 
parametrization of the hadronic matrix element for J/^lJ to vector meson given in Eq. ([3]) is less studied 
before. Similar matrix element for a transition of a vector to another vector meson which is induced by 
the electromagnetic current was investigated by Kagan in Ref. [itI . 
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B. The transition form factors in the QCD sum rules 

(*) — 

In this subsection, we calculate transition form factors of J — > Df^^ s-^ iii the QCD sum rules. Here 
we present the formulations for J /ip ^ oi*^ transition explicitly, while the expressions for J/ip D^*^~ 
can be obtained by simple replacements of D^*^ D^*^~ and s quark to d quark. 



The matrix element for J /ip — !■ 



Following the standard procedure of the QCD sum rules [9[, we write the three-point correlation 
function for J /ip to Dj as 

n^, = [ d''xdSe-^P^ y+'P'--{0\j^^ {x)j^{0)j;^/^{y)\0), (4) 



where the current ju^^{y) = c{y)'^i,c{y) represents the J/if) channel; ^/^(O) = 57^(1 — 75)0 is the weak 
current and (x) = c{x)i^^s{x) corresponds to the D~ channel. In terms of the following definitions, 

(0|c7,c| = m^Uel, (0|cf75s|D.) = I^fll^^ (5) 

rric + rris 

we can insert a complete set of hadronic states with the quantum numbers the same as J /if: and Dj to 
achieve the hadronic representation of the correlator ^ 



fD,ml)^{Ds\j/,\J/i^)m^f^ 



Uuu = 7 — ^— ^ — /"l^' — ' '^L Z^'"'' ^ — - + contributions from higher states. (6) 

Obviously the concerned lowest hadronic states are J/ip and Dg, while the terms with "higher states" 
represent contributions coming from higher excited states and continuum. Using the double dispersion 
relation, the contributions of excited states and continuum can be expressed as 

contributions from higher states = / / dsids2 , ^^'^^ nV/ ^''^ ^ + subtraction terms, (7) 

where X]i2 denotes the integration region in the {si,S2) plane. is the spectral density at the hadron 
level. The subtraction terms are polynomials of either pi or which should disappear after performing 
the double Borel transformation B]^j2Bjyj2, with 



(-p2)(n+l) / ^ y 

nl \ dpf 

-pf/n=M^ 



^M?= (312 ■ (8) 



On the other side, we calculate the correlation function at the quark level by using the OPE as 

^liiy = -foe^ual3PlP2 " + f2P2f,P2y + hP2fj,Ply + fm,,P2u + hOl^u), (9) 
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where each coefficient contains contributions from both perturbative part and the non-perturbative part 
whose effects manifest in several typical condensates, 

/, = /r*I + frim) + f?''{GG) + ff^qGq) + (10) 

with /f^^*, ff^, fi^'^i ■■■ denoting the contributions to the correlation functions from dimension 0, 3, 
4, 5,... operators. By the quark-hadron duality, one may match the two different representations of the 
correlation function and perform the double Borel transformation on variables pi and p2, then we get the 
sum rules for the form factors 

V{q') = - + "^-^^^"-^ \'^''^\-^V^'h<s/^'HllMlB /o, (11) 

(m^ + mDsjm^UjDsmjj^ 
^^(^,) ^ _ (m + mjirn, + rnnJ ^^l/Mi^rnlJMl^2^2j^^^ ^ ^^^^ (^3^ 

2. The matrix element for J/ip D*^ 

The three-point correlation function of J/ip to D*~ is 

n^.p = i'J d^xdSe-'P^-y+'^^-^iS)\j^*^ (x)i^(0)i,^/'^(2/)|0), (15) 

where the current j^" {x) = c{x)^ps{x) denotes the D*^ channel, and ju^^ (y), jfi{0) are defined as in the 
above subsection. One the one hand, inserting the hadron states, the correlation function is written as 

n - "^^^-^^^-^p''^^-^'^'^'-^/^)"^^/^^^/^^^' + f fds, + subtraction terms (16) 

(m^^ - pj)iml, - Pi) + y y '^"^^'^ {s^ - pl){s2 - pi) + terms.(lb) 

On the other hand, the correlation function at the quark level is formulated as 

n^i^p = 'iFl^iiua(3PlP2Plp + 'iF2ef,ual3PlP2P2p + iF3e^pai3PlP2Plu + iF4^fipal3PlP2P2u 

+ iF^^vpcci3PiP2Pip. + 'i'Fee„pai3PiP2P2p + F-jg^^pip + Fsg^pPi^ + Fgg^ppi^ + FiQg^yP2p 

+ FngppP2u + Fi2gupP2p + Fi^iPi^pi^pip + Fi4P2pP2^Plp + Fi5PipP2^Plp + Fi(iP2pPl^Plp 

+ Fl7P2pP2uP2p + FisPlpPl^P2p + FigP2^Pl^P2p + F20P1 ^,P2uPl p, (17) 

where each coefficient Fi includes contributions from both perturbative and nonperturbative parts, and 
is written explicitly as 

Fi = *I + + i^«^(GG) + Ff\qGq) + .... (18) 
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Again, equating the correlation functions calculated in these two frameworks and performing the Borel 
transformations on both sides, we derive the form factors ol J/ip ^ D*~ as 

1 



%{q' 

Mq^ 
Ve{q^ 



ml 



iLe</Ml^<tlMiMlMlB{F, + F, 



1 



M(MiB{Fi + F( 

2 /A/r2 



6 



e™yM^^-i,j/M|^2^2^^^^ ^ ^^^^^ 



(20) 

(21) 

(22) 
(23) 
(24) 

(25) 



2m^/^mz)*/Dj 

^ e<"'"^e'^"-t''''^ MIm'MFi, - F,,)q' + (Fm + F,,){ml. - m/)], (26) 

1 .™yMf -i*/M| 2,.2 



e -^/-^i e -D- - MtM^BFii, 



1 



MfMiBFj 



(27) 
(28) 



C. The Wilson coefficients for correlation function Ylni, 



In this subsection we calculate the Wilson coefficients which are defined above. To guarantee sufficient 
theoretical accuracy, the correlation functions are required to be expanded up to dimension-5 operators, 
namely quark-gluon mixing condensate. The dimension-6 operators, such as the four quark condensates, 
are small and further suppressed by O(a^), so can be safely neglected in our calculations. 

The diagrams which depict the contributions from the perturbative part and nonperturbative con- 
densates are shown in Fig. [TJ The first diagram results in the Wilson coefficient of the unit operator; the 
second diagram is relevant to the contribution of quark condensate, where the heavy-quark condensate 
is neglected. The Wilson coefficient of the two-gluon condensate operator is obtained from Fig. [U^c-h). 
The last two diagrams Fig. HJi-j) stand for the contribution of quark-gluon mixing condensate. In this 
work, all of the Wilson coefficients are calculated at the lowest order in the running coupling constant of 
strong interaction. 
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(d) (e) (/) 




{g) ' ih) (^) 




U) 



FIG. 1: Graphs for the Wilson coefhcients in the operator product expansion of the correlation function, (a) is for 
the contribution of unit operator; (b) for the two-quark condensate; (c-h) describe the contributions from giuon 
condensate, (i-j) is for the quark-gluon mixing condensate. 



1. Perturbative contributions to Wilson coefficients for Hfj^u 

The perturbative contribution to the three-point correlation function li^u shown in Fig. [T] (a) is 
included in the following amplitude 

= I T^ti (-l)'I^[>ir^^^5--4 7.(1 - 75)-— 4 ]> (29) 

where rrtq denotes the mass of the light quark in the D meson, and the factor "3" is due to the color 
loop. Using the dispersion relation, C"^^^ is written as 
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The integration region is determined by the following condition 

_ ^ ^ 2gi(s2 + - ml) - si{si + S2 - q^) ^ ^ ^^^^ 

where A(a, b, c) = + b'^ + — 2ab — 2ac — 26c. The standard way to calculate the spectral function 
Piiu{si-,S2-,q^) is described below 0]: Firstly, it is essential to calculate the double discontinuity of the 
amplitude, which can be realized by putting all the internal quark lines of Fig. [T] (a) on their mass-shell 
and substituting the denominators of the quark propagators by the 5 functions based on the Cutkosky's 
cutting rule, 

-5 \ > -l-nidik^ -Tn^). (32) 

— + le 

Then, the spectral function can be easily achieved. Finally, we get the expression of the spectral function 
in the following form 

5{k^ - ml)5[{p2 + kf - ml]6[{pi - kf - ml]. (33) 

After tedious calculations, one finally obtains the perturbative contribution to the correlation function, 
which can be decomposed as the sum of various terms according to different Lorentz structures, namely, 

f^iZ^ = -Po'''^e^lya|3PtP2 - ^{f>T''^Pl^lPlu + pT'^*P2f^P2u + pT^P2^,PIu + pT^Pl^lP2u + pT'^^a^iu)- (34) 

The expressions for the pf^^* are a bit more tedious, so that we will display their explicit forms in 
Appendix Rl 

2. The quark condensate contribution 

Now we turn to calculating the Wilson coefficient of the quark condensate operator, which are shown 
in Fig. [11(b). One can easily find that it does not contribute to the correlation function after performing 
the double Borel transformation on both variables pf and since the propagator of this diagram 

■j-^ JTT^ 2T only depends on variable p^- other words, the Wilson coefficient of dimension-3 two 

quark condensate turns to be zero in the leading order of heavy quark mass expansion after carrying out 
the double Borel transformation. As can be seen, vanishing of the contributions from quark condensate is 
independent on the structures of the effective vertices, therefore, it also does not contribute to the decays 
of J/ip into a vector meson for the same reason. Below, we do not need to investigate the contributions 
of quark condensate to J/tp ^ D* based on this argument. 
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3. The contribution from gluon condensate 

The diagrams which determine the Wilson coefficient of the gluon condensate are shown in Fig. [1] 
(c-h). The standard way is using the so-called fixed -point gauge technique. The gauge fixing condition 
is 



x^^^ = 0, 



(35) 



where A" is the gluon field. In the momentum space, A"j^{k) is transformed to the gauge invariant field 



strength as 



(2^)*G^^(0)— 5^(^) + ... 



dkn 



(36) 



Indeed, the loop integral 



.Mn(a,^> c) 



d^k 



(37) 



(27r)4 [A;2 - m2]«[(pi + A;)2 - mf]^[(p2 + k)^ - mlY' 

which is encountered in the work, is not easy to be performed by the Feynman parameter method. On e 
alternative way to calculate this kind of integrals has been extensively discussed in the Ref. 



3, 



i8i,Q,y], 



where the authors suggested to work in the Euclidean space-time and employ the Schwinger representation 
for propagators. Instead, in our work, we follow the method employed in Ref. 2l|, [22], namely, directly 
calculate the imaginary part of the integrals in terms of the Cutkosky's rule. 

With the help of the Mathematical package "FeynCalc" , we finally get the contributions of Fig. [1] 
(c-h) at the price of some long and tedious derivations and time-consuming computer computations. The 
contributions of the gluon-condensates from various sources cancel each other completely after carrying 
out the double Borel transformation to the variables Pi and p^- Therefore the diagrams involving the 
gluon condensate do not contribute to the transition of vector meson J to a pseudoscalar D meson. 
This argument also applies to the transition of a pseudoscalar meson to a vector which was discussed in 



Ref. 



21 



22( 1 ■ since topologies of the Feynman diagrams which result in the Wilson coefficient of the gluon 



condensate are the same. As analyzed later, it is also true for the transition of J to a vector meson. 
However, we find that the flavor-changing neutral current process can receive non-zero contributions from 
the gluon condensate. It should be noted that the null contributions of g. 
for the weak transition c — > s{d) are different from that obtained in Ref. 



uon condensates to sum rules 

nn 



m Il9|,l20|, where the method 



they adopted does not allow for the substraction of continuum contributions. 



10 



4- The quark- gluon mixing condensate contribution 

Finally, we go on calculating the Wilson coefficients of the dimension-5 operator (qGq). Only two 
diagrams shown in Fig. [1] (i-j) are involved. Concentrating on these two diagrams, we find that they do 
not contribute to the correlation function, due to the same reason as that for the null contribution from 
quark condensate, namely, only the variable appears in the propagators, the amplitude will vanish due 
to the double Borel transformation. 

As mentioned at the beginning of this section, we do not consider the four quark condensate, hence 
only the perturbative part which corresponds to Fig.l (a), offers a non-zero contribution to the correlation 
function. 



D. The Wilson coefficients for the operators contributing to the correlation function n^j.^ 

After above lengthy discussions, a computation of the correlation function H^up which determines the 
transition amplitude of J ftp to a vector meson is straightforward. Repeating the previous calculations 
but replacing the vertex for the pseudoscalar meson to that for a vector meson, one can obtain the 
expressions of the Wilson coefficients for all the concerned operators. 

1. The calculations of the perturbative contribution to Tlf^iyp 
The Wilson coefficient of the perturbative part corresponding to Fig. [1] (a) is 

We rewrite it in the form of dispersion integrals for the sake of connecting it to the hadronic spectral 
density based on the assumption of the quark-hadron duality, as 

The integration region is the same as that for the C^^^*, which is presented in Eq. (j3ip . Setting all the 
internal quark lines on their mass shells, we derive the spectral function p^fj^ as 

nert ■ /pert a , ■ ipert ad , ■ /pert a , ■ /pert a 

P]iup = ^P l <^pua0PlP2Plp + 1'P 2 epua0PlP2P2p + iP3 ^tJipa0PlP2Plv + ^P A ^ppa0PlP2P2u 

I • /pert a , • /pert a , /pert , /pert , /pert 

+IP5 (^upa0PlP2PltM + I'P 6 <^upa0PlP2P2p + P7 df^uPlp + Ps 9ppPlu + P 9 9upPl^ 

, /pert , /pert , /pert , /pert , /pert , /pert 

+PlO 9f^uP2p + P II 9f,pP2u + Pl2 9vpP2^, + Pvi PlpPluPlp + P li P2^P2,yPlp + P 15 PlpPluPlp 

I /pert , /pert , /pert , /pert , /pert 

+PW P2pPluPlp + P 17 P2pP2uP2p + P Ifi PlpPlvP2p + Pl9 P2f^PluP2p + P 20 PlpP2uPlp- (40) 
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Only /j'f * {i = 1, 4, 5, 6, 7, 9, 11, 12, 14, 15) are related to the form factors Vi, Vz, V3, V4, V5, Ve, ^1, ^, 
^3, ^4, and we display their expressions in Appendix [Bl 

^. The contribution of gluon condensate to II^i^p 

Similar to the derivation made above, we easily obtain the Wilson coefficient of the gluon condensate 
which may contribute to the correlation function 11^,^^. Then we rewrite the Wilson coefficient in the 
form of dispersion integrals: 

= / / d^.d,. /"-"'^;-'""'' . (41) 

^ J J {si -pi)[s2-pi) 

where the integral region is the same as that for the perturbative part. 

The Lorentz structures corresponding to pl (i = 1, 4, 5, 6, 7, 9, 11, 12, 14, 15) are 

GG ■ iGG \ , . ,GG X , . /GG a P , ■ iGG a P , /GG 

Pfiup = ^Pl WAPl+«P4 WAP2+^/'5 ^l^uaPPlPlPlu + l-P 6 f^tipafiPlPlPlp + P 7 9puPlp 

GG GG GG GG GG 

+P'9 9upPlp + P'll gppP2y + P'l2 9upP2p + P'l4 9ppP2y + P'l5 9upP2p + ■■■■ (42) 

After some long and tedious calculations, we find that all of the above coefficients p'f *^ are zero. This is 
completely the same as for the 11^,^ case. Therefore, only the perturbative part survives after performing 
the double Borel transformation on the two variables p\ and p'2 at the leading order of the heavy quark 
mass expansion and QCD running coupling constant expansion for the three-point function Ilfj,up- 

III. NUMERICAL RESULTS OF FORM FACTORS IN QCD SUM RULES 

ate form factors numerically. Firstly, we collect the input parameters used 



Now we are able to ca' 
in this work as below 23, 



cu 



24 



25| 



mc(mc) = 1.275 ± 0.015GeV, m^(lGeV) = 142MeV, 
m„(lGeV) = 2.8MeV, md(lGeV) = 6.8MeV, 

a^(lGeV) = 0.517, nij/^ = 3.097GeV, 

nir)- = 1.869GeV, m^- = 1.968GeV, 

(43) 

niD*- = 2.010GeV, m^.- = 2.112GeV, 

fj/^ = 337j:|2MeV, fo- = 166l?oMeV, 

= 189l?oMeV, /b.- = 240lJ°MeV, 

/^,_ = 262t?2MeV. 

All the QCD parameters are adopted at the renormalization scale around 1 GeV. It should be pointed out 
that the mass of charm quark used in this work is determined form the charmonium spectrum in Ref. 



23|. 
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As for the decay constants of charmed mesons, on the one hand, there is a flood of papers on the theoretical 



27 



28|, 



2S 



3C 



31 



32 



33 



34 



35, ea,^ 



on 



investigation of leptonic decay constants of and Dg j26l . 
the other hand, the measurements of decay constants of pseudoscalar and Ds mesons have recently 
been improved by the CLEO and BaBar collaborations [38l . |39|. Moreover, the CLEO collaboration 



reported their work on the value of ratio fr,+ /fD+ using the measurement of Df I'^v channel and 



41 



42l |. However, the decay constants of D*^ and D* mesons have 



obtained /„+ = 274 ± 13 ± 7MeV 

not been directly measured in experiments so far. The only available results on //)*+ and f]j*o from 
the Lattice QCD calculations 



29|, Q, Q determine fo* 



272 ± leiJnMeV that is smaller than the 



value of decay constant for D'^ measured by the CLEO collaboration 



41 



42| . To reduce the theoretical 



uncertainties in the three-point sum rules of the weak transition form factors, due to quarks masses, 



threshold parameters and Coulomb-like corrections of J /ip effectively [43|, we use the decay constants 



and /^(*)- calculated from the two-point QCD sum rules in leading order of a^, the same as that in the 



three-point sum rules. The explicit calculations of the decay constants, in the framework of QCD sum 



rules, for both J/il) and L*^*] are displayed in Appendix O Our results indicae tha t U 



1.1, 



42|. 



which are in good agreement with that from lattice simulation [40|] and experiments 41 

For the threshold parameters and one should determine them by demanding the QCD sum 
rules results to be relatively stable in allowed regions for Aff and Aff , the values of which should be 
around the mass square of the corresponding first excited states. As for the heavy-light mesons, the 
standard value of the threshold in the X channel would be = {nix + ^x)'^, where Ax is about 0.6 



GeV 



44 



45 



46 



471, 



and we simply take it as (0.6 it 0.1) GeV for the error estimate in the numerical 



45|, 



46 



we select 



analysis. When it comes to the heavy quarkonium, following the method in Ref. 
the effective threshold parameter to ensure the appearance of the pleasant platform and also around the 
mass square of ^'(25'). In this way, the contributions from both the excited states including i/j{2S) and 
the continuum states are contained in the spectral function. 



A. The numerical results of the form factors 

1. Evaluation of the form factors for the J/tp D~ 

With all the parameters listed above, we can obtain the numerical values of the form factors. The 
form factors should not depend on the Borel masses Mi and M2 in a complete theory. However, as we 
truncate the operator product expansion up to dimension-5 and keep the perturbative expansion in to 
leading order, an obvious dependence of the form factors on these two Borel parameters would emerge. 
Therefore, one should look for a region where the results only mildly vary with respect to the Borel 
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masses, so that the truncation is reasonable and acceptable. 

With a careful analysis, Si = 13.7 GeV^ and S2 = 6.1 GeV^ are chosen for the form factor V calculation. 
We require the contributions from the higher states to be less than 30 % and the value of V does not 
vary drastically within the selected region for the Borel masses. As commonly understood, the Borel 
parameters and M| should not be too large in order to insure that the contributions from the higher 
excited states and continuum are not too significant. On the other hand, the Borel masses also could 
not be too small for the sake of validity of OPE in the deep Euclidean region, since the contributions of 
higher dimension operators pertain to the higher orders in -Yrii = 1,2). Different from that adopted in 

n n 

previous literature [111 . |15I | where the ratio of Mi and M2 was fixed, in the calculation of form factors, 
we let Ml and M2 vary independently as suggested by the authors of Ref. Q, [3]. In this way, we 
indeed find a Borel platform Mf G [6.0, 10.0]GeV^M| G [1.0, 2.0]GeV^ plotted in Fig. [21 which satisfy 
the conditions discussed above. One can directly read from this figure that V{q^ = 0) is OASto^^l, whose 
uncertainties originate from the variation of the Borel parameters. 

Following the same procedure, we also obtain numerical results for the other three form factors Aq, Ai 
and A2 within the chosen Borel window as shown in Fig. [2j The numerical results of form factors V, Aq, Ai 
and A2 at zero momentum transfer are then 

V{0) = OMt'^.Z Ao{0) = 0.27l[i:|]2, Ai(0) = 0.27t[J;°3, ^[^(O) = 0.34;°;°^. (44) 

It needs to be emphasized that the form factors ^i(g^), ^2(9^) and AqIq'^) should satisfy the relation 
(m^ + m^-)Ai(O) + {m-^ — mj^-)A2{0) = 2m^AQ{0) to ensure disappearance of the divergence at the 
pole = 0. The theoretical uncertainties in the form factors (j44p originate from the Borel masses 
and Ml . They are at the level of 15%, which implies stable results from the QCD sum rules approach. 

Indeed there are some extra errors originating from the values of Si and 53 which correspond to the 
threshold of the higher excited resonances and continuum states for the J ftp and D channels respectively. 
In the QCD sum rules approach, the values of threshold parameter is usually in the vicinity of mass square 
of the first physical excited state, therefore, we do not investigate the dependence of form factors on the 
threshold parameter in this work as that in Ref. [l^, [3] 1 where a larger threshold value of charmonium 
is adopted. This uncertainty would cause errors in the resultant form factors. Besides, the fluctuations 
of the charm quark mass can also result in the uncertainties of the form factors, which is evaluated to 
be at the level of 6 — 8%. Moreover, the input parameters such as the decay constants of D meson and 
J/ip can also bring on additional uncertainties. Combing the errors from various parameters discussed 
above, the uncertainties on the form factors can be estimated within 20 to 30%, expected by the general 
understanding of the theoretical framework. 
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FIG. 2: Dependence of form factors V,Aq,Ai and A2 at = responsible for the decay of J/i/j D on the 
Borcl masses. 



Next, we can further investigate the dependence of the form factors V, Aq, Ai and A2. The physical 
region of for J/ip D^l'^ui is < < {ttij/^ — m^)-)^ ~ l.SGeV^. However, with the QCD sum 
rules, we could not obtain the form factors in the whole physical region, since the additional singularities 
- so called "non-Landau-type" singularities emerge, which had been extensively discussed in Ref.[ll|. To 



avoid this kind of singularity, we restrict our calculations in the range 

of g2 g [0,0.47]GeV2. We show 

the dependence of the form factors V, Aq, Ai and A2 in Fig. [3l 
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FIG. 3: dependence of form factors 1/, Ao,j4i and A2 for J/V-J 
non-Landau-type singularities. 



0.4 



D within the kinematical region without 



In addition, for the convenience of applications to phenomenology, one can parameterize the above 



form factors in a three-parameter form 



F^il') = 2, ^/^L 4/ 4 . (45) 

where Fi denotes the form factors V, Aq, Ai and A2, and Oj and bi are the parameters to be fixed. Using 
the QCD sum rules Fi{q^) with q'^ restricted within a certain kinematic region, we can fix the parameters 
Cj, bi in the expression. This double-pole expression for form factors can be generalized to the whole 
kinematic region. Finally, our results for the parameters Oj, bi are given as 



aAi 



0.931°:?^, bA, = 0.46t0:g?, aA, = lAltlH ^a, = ^■^2tlfi- (46) 



For the other form factors which are discussed in the following subsections, we will adopt the same 
procedure to obtain the form factors in the whole kinematic region. 



2. J lijj — > form factors 

Now, we move on to the computations of the form factors for the transition J/ip — > D~ , which is 
quite similar to that for J/ip , only with d quark in D~ being replaced by s. It is also noted 

that the threshold parameter = 6.6GeV^ for the Ds channel and the Borel window are shifted slightly 
compared with that of J/V' — > D~ . Since the figures are very similar to the case for J/^lJ — > D~ , we just 
omit them. The obtained form factors for J/tp D~ at q^ = are 

y(0) = imtl-'ol AoiO) = 0.37t0;0i, ^^(O) = 0.381°:°?, ^2(0) = 0.35t[|;08. (47) 
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FIG. 4: dependence of form factors V^Aq^Ai and A2 for J/ip within the kinematical region without 

non-Landau-typc singularities. 

And the parameters Oj and hi defined above for the dependence formula (with the replacement of 
D ^ Ds) are fixed as 

ay = 1.86l°;26, by = Om+'^^f^, a^, = 2.12l0;04, = l.SOll];^^, 

UA, = l.lSt^om, bA, = O^rf'oZ aA, = lAlt'oil, bA, = O.SSt^ofv (48) 

Again, the relation {m^ + m£)JAi{0) + {m^ — m£)JA2{0) = 2m^^o(0) is well respected, which guarantees 
that the hadronic matrix element responsible for the J/ip ^ Dj is free of divergence due to the pole at 
= 0. We show the dependence of the form factors on in Fig HI 

3. J/ip — > D*^ form factors 

The evaluation of the form factors responsible for J /tp ^ D*~ is performed following the standard 
procedure, with appropriate Borel windows obtained. The threshold value for D*~ channel takes S2 = 
6.8GeV^ in our numerical analysis. The form factors at zero momentum transfer q"^ = are collected 
below as 



Alio) 


— U.4U_o.oi) 


i2(0) 


- n 44+0-10 


i3(0) 


= 0.861°:°?, 


^4(0) 


_ nn, +0.06 

— u.yi_o.o4> 






V2(0) 


— U.UO_Q Q4, 


^^3(0) = 


_ n 99+0.03 


V4(0) 


— U.ZD_QQ5, 


^5(0) 


_ 1 07+0.08 

— -'^•■-''-0.03' 


V6(0) 


- n 87+0-05 
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From the above results, we find that the form factors obtained in the QCD sum rules respect the relations 
^i(O) = A2{0) and V3(0) = V4(0), which are essential to assure that the hadronic matrix element of 
J/if) ^ D*~ is free of divergence at q'^ = 0. 

Different from that discussed for the J/ip ^ ^ds case, not all the form factors which appear in the 
hadronic matrix element for J/ip D*~ , are suitably parameterized in the form of eq. ()45p with the 
thee-parameter approximation. To be more specific, the dependence of the form factors Ai{q^) and 



^2(9 ) are written in the following form 5ll. |52| 



where Fi represents Ai and A2; while the other eight form factors are written in the three-parameter 
form, 

G^i<^') = 1 4/ 4 . (51) 

where Gi can be ^3, A4 and Vi{i = 1 — 6). We then extend the form factors to the whole physical region 
< q^ < {mji^ — rriu*-)^ ~ 1.2GeV^, by fitting the parameters as 



Oil 


_ 1 77+0.04 
— '-0.01' 


«A2 = 


- 1 qt;+0.17 

- i.yO_o.25i 


"A3 


- 2 qS+°-^S 

— z.yo_Q 




_ 9 47+0.54 
— ^•^'-0.27' 




9 70+0. 05 

— Z./0_g Q3, 




1 70+0.27 




_ 1 qfi+0.03 




_ n 08+0.07 

— u.y5_Q Qg, 




_ 9 n+0.04 

— Z.ii_QQ4, 




"-'•^J^-0.02; 




- 1 q2+°-° 

— i.JZ_QQ3, 




— 1 C7+0.12 


«y4 


- ^■yo_o.23' 




1 07+1-03 
J^-y ' -0.34' 




- 1 q2+°-° 




- 1 ns+°-° 

— i.U0_Q2i, 


«V6 


— Z.UU_o.o2; 




J--UO-0.03- 











(52) 



4-. The form factors for J/ip s- D* 

The computation on the amplitude oi J/ip ^ D*~ is almost the same as that for J/'i/; — > Z?*~, only 
d quark in D*~ being replaced by s quark, with the difference resulting in a different Borel platform. 
Besides, the threshold parameter for the D*~ channel is set as S2 = 7.4GeV^ in the calculations. The 
q dependence of the form factors falling into the region of q^ e [0,0.37]GeV2 is plotted in Fig. El As 
mentioned before, the form factors Ai and A2 can be parameterized in the form of Eq. (j50p . while the 
other form factors can be fit in the usual three-parameter form in Eq. ()5ip . The parameters Oj and bi can 
be determined by reproducing the numbers obtained from the QCD sum rules for the kinematic region 
q^ G [0,0.37]GeV^ and then we generalize the results to the whole physical region q^ G [0, 0.97]GeV^. 
The values of these parameters together with the form factors at = are collected for convenience as. 
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FIG. 5: dependence of form factors Ai, A2, A3, ^4, Vi, V2, V3, V4, V5 and Vg for J/i! D* within the kinematical 
region without non-Landau- type singularities. 



and 



„ _ 1 Q9+0.0 _ _ -1 oc+O.Ol ^ _ q r,7+0.12 . _ 1 QR+0-80 

"Ai ~ -"^-^^-o.os' "^2 ~ ^•"'-'-0.21' "^3 ~ •-'•^'-0.01' "A3 ~ -"--^"-o.ie' 

"^^4 = 3.08j:Q;o2, = 2.08j:o;26, = 2.05to;o2, Vi = O-QO^aoe' 

ay, = 2.53lHi by^ = Omtlil a^^ = 2.0<1^, by^ = 2.1<o^ 

„, _ q 09-^0.37 i_ _ 1 ye -HI. 58 _ 1 09-^0.22 7_ _nQ-|-H0.40 

"Vi ~ "^-"^^-O.aS' V4 ~ -'^•'"-0.67' "Vs ~ -"^-^^-O.OS' V5 ~ '-'•"-'--o.os' 

ay^ = 2mtlil by^ = 0.8lt°:i, (53) 



ii(0) = 0.53t|];0?, i2(0) = 0.53tHf , A,{0) = OMf'o'ol A,{0) = 0.9lt|];[]?, 
Vi(O) = 0.5<oi, 1/2(0) = 0.691°:°^, ^3(0) = 0.2<o?, UO) = 0.261^^, 
F5(0) = imt'oll MO) = l.l<o?. (54) 

In the same way, the relations Ai(0) = ^2(0) and V3(0) = ^4(0) are well satisfied. 



IV. DECAY RATES FOR SEMI-LEPTONIC WEAK DECAYS OF J/V^ 



With the form factors derived above, we can perform calculations on partial widths of the semi-leptonic 
decays of J/ip. The relevant CKM parameters are directly taken from the particle data book 24i |: 



iKrfl = 0.2271, 



0.973. 



(55) 
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FIG. 6: dependence of form factors Ai, A2, A3, A4, Vi, V2, V3, V4, V5 and Vg for J/i^ ^ D* within the kinematical 
region without non-Landau- type singularities. 



as 



For the semi-leptonic decays J/ip — > D^^^^l'^i'i {I = e,/i), the differential partial decay rate is written 
dT 



1 1 



3 (27r)3 32m3 7«„„ 



|M, „(.)-,, I du, 



(56) 



dg2 

where u = {pi+ +Pui)'^', Pi+ and p^i are the momenta of and vi respectively; |Mp is square of the 

, (*)— 

transition amplitude after integrating over the angle between the r and D^^^y The upper and lower 
bounds for u are given as 



mf, — m 



D 



(•) 

(d,s) 



my 



(57) 



The transition amplitude M for J/^p — > D^*^^-^l'^ui reads as 



M, 



Gp 



Finally we get the branching ratios of the semi-leptonic decays as 



(58) 



BR{J/^ - 


D- 


e^Ue) -- 


= 7.3^2:^ X 10- 


12 


BR(J/V' - 


D- 






= 7.lt|2 X 10- 


12 


BR{J/^ - 


-D- 


e+Ue) -- 


= 1.8l°i X 10- 


10 


BR( J/V' - 








= i-7iL' X 10- 


10 


BR(J/V' - 


D* 


-e+i/e) 


= 3.7l|;? X 10 


-11 


BR(J/V' - 


D* 






= 3.6^1:? X 10 


-11 


BR(J/?/; - 


-* D* 

^ s 


"e+i/g) 


= f>.Qt\l X 10 


-10 


BR(J/V' - 


D* 






= f>At\l X 10 


-10 



(59) 



where we have combined various uncertainties in the form factors discussed in last section to determine 
the final error tolerance in our theoretical calculations. Our predictions are much below the present 
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experimental upper bounds BR{J/ip D~e^i^e + c.c.) < 4.9 x 10~^, BR(J/V' D~e^Ve + c.c.) < 
1.2 X IQ-^ 

A few remarks are presented in order. Firstly, the sum of the branching fractions of semi-leptonic 
decays of J ftp whose final state includes D~ , D*~ , e and fi and their charge conjugate channels can 
reach as large as 3.1 x 10"^, which is expected to be marginally observed at BESIII. Secondly, it is 
worthwhile to point out that the decay rates for the dominant semi-leptonic weak decays of J/V' obtained 
in Ref. [l| were about 7 x 10~^, which are two times greater than that calculated in this work. This 
discrepancy can attribute to the heavy quark spin symmetry and the non-recoil approximation used in 
Ref. [l|], also to the different methods used to estimate the non-perturbative form factors ^. Thirdly, the 
ratio of i?i = ^^^f(f~^^^^ ^, ^''} ~ 3.1 is about 2 times larger than the value calculated in [l|, where the 

BR( J/i/i^Dg e+i/e) " 

assumption of heavy quark spin symmetry and non-recoil approximation were adopted. Forth, the ratios 
- l^iJ^Z'^nS^]) R, ^ lllf^Z'S-'Sli should be equal to ^ 18.4 under the SU(3) flavor 
symmetry limit. Our numerical calculations show that R2 — 24.7 and ~ 15.1, which implies large 
effect of SU(3) symmetry breaking. 

V. DISCUSSIONS AND CONCLUSIONS 

The charmonium J/ip meson can decay via the strong and electromagnetic interactions, thus weak 
decays of J/ij: should be very rare unless there is new physics beyond the standard model to make a 
substantial contribution. If such weak decays can be measured by the future experiments with sizable 
branching ratios, it would be a clear signal for new physics. 

To make the new physics signal clearly distinguishable from the standard model, a careful study 
of weak decays in SM is needed. In this work, we calculated the form factors of weak transitions of 
J/il) ^ D{^(is) terms of the QCD sum rules. With the form factors, we estimate the branching ratios of 
the semileptonic weak decays of J /'ip and find that the sum of the branching ratios corresponding to the 
dominant modes is about 3.1 x 10^^ which may be marginally measured by BES III. The QCD sum rules 
approach possesses uncontrollable errors as large as 20 to 30%, was confirmed by our numerical results. 
Moreover, due to a Coulomb- type correction in heavy quarkonium decay (or i?c)j which may manifest as 

OOP diagrams), the spectral function needs 



the ladder structure in the loop-triangle (as a part of multr 
to be multiplied by a finite renormalization factor ^. 



191, 



541 ]. This would bring up another kind of 



uncertainty. It is expected that this kind of correction can give birth to the twice multiplication of the 



1 In Ref. d, the ISGW model [53| was employed to compute the single form factor 7^12, while we adopt the QCD sum rules 
to calculate the form factors in this work. 
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form factors at the maximal momentum transfer. However, in this work, we calculate both the three-point 
QCD sum rules of the weak transition form factors and the two-point sum rules for the decay constant 
of J/ip to the same order of a^- Then, it is expected that most uncertainties due to the Coulomb-like 
corrections are canceled in our calculations; therefore, the Coulomb-like corrections for the J /ip channel 
are not included in our calculations. As for the heavy-light mesons, there are no corrections in the power 
of inverse velocity for it, since the light quark moves relativistically. Therefore, one should explore the 
sum rules for both three-point and two-point correlation functions up to next-to-leading order in strong 
coupling constant so that Coulomb-like corrections to the heavy-light mesons can be canceled effectively. 
Moreover, an explicit calculation of Coulomb-like corrections to the heavy- light vertex in triangle diagram 
is still not available now, which can be left for further considerations. 

One can trust the numerical results to a certain accuracy, at least the order of magnitude is reliable. 
With these form factors we may continue to estimate the rates of non-leptonic weak decays of J /ip as 



55|. 



long as the factorization theorem is proved. That would be the contents of our next work 

The branching ratios of semi-leptonic weak decays of J /ip are very small in SM, even though their 
strong decay modes are OZI-suppressed. Our numerical results indicate that even with a large database 
which will be collected by the BES III, such weak decay modes may only be marginally observed. There- 
fore we lay our expectation on our BES III colleagues and hope them to provide sufficiently large database 
to make this challengeable field more fruitful. 
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APPENDIX A: THE EXPLICIT FORMS OF WILSON COEFFICIENTS FOR H^^ 

In this appendix, we would like to show the explicit expressions of Wilson coefficients appearing in 
Eq. ljimH]) after the Borel transformation. As mentioned before, only perturbative part contributes to 
the correlation function for J/ ip decays to D^^ at the leading order of heavy quark mass expansion and 
the Og expansion series. Namely, Eq. (jlO|) can be written as 

/, = /r*I + 0(a,) + 0(l/m,), (Al) 



22 



pert 



with h being the heavy charm quark here. The /f^^* can be related to p^^^^ defined in Eq. (j34p by 

'(mc+rriq)' 

or 



dS2 I dsi- ^ 



pert 



dS2 



1 g~sl/A/2^ 



Mr 



-S2/MI pert 



Pi 



(si,S2,g^)- 



The lowest bound of si, i.e., sf can be determined by the Eq. ([HT|) as 
1 



2mg 



(2m + S2 + 9 + m + S2q - rn {s2 + q ) 



+ y?7^c ~ 2(m2 + S2)m1 + (m^ — S2)^ym^ — 2(m^ + q^)m1 + (m^ — (72)2 
according to the Landau equation 



The obvious forms of 



pert 



0,2,4,5) are 



Po'"'*(si,S2,g^) 
f>2''\si,S2,q^) 



■^-y^^[mcA + (mc - mq)si{2ml - Im^ - si + S2 + g^)], 
3si 



2^2A5/2 



(rUc - mq)si 



6m^ — 6(2m„ + si — S2)m^ + 6m + (si — S2) 



+6mg(si - S2 



pert/ 2\ 

Pa {si,S2,q ) 



(A2) 



(A3) 



(A4) 



pert I 2 



+ mc[2(?nc — mq){2mc + m^) — si + S2]A 

+q'^[2(mc - mg)si(3m^ - 3mg - si + 2^2) + rricX + {rric - mg)sig^]|, 

4^2*^5/2 { [~"^cA^ + {rric - mg){2s2ml + si(2mg + si - S2)]A 
+2(mc - mg)si[3(si + S2)"^c 

-2(3(si + S2)m^ + {si - S2)(si + 2s2))m^ + (si - ^2)^52 + 3mg(si + S2) 
+2m^(si - S2)(si + 2s2)] + (mc - mg)g2[2si(s| + (-2m^ + 2m^ + si)s2 
+{ml - ml){-^ml + 3m2 + 4si)) - (2m2 + si)A - 4si(m2 - m^ + S2)g^]|, 

A(mqSi + mcS2 — rricq^) — 2(mc — mq)[Am^ 



87r2A3/2 



+ (m^ - mg)si(77T,^ - m^ - si + S2) + si(m^ - m^ + S2)g^]|. 
In this appendix, we adopt the notion A = A(si, ,52, q'^) for the convenience of writing. 



(A5) 



APPENDIX B: THE EXPRESSIONS OF WILSON COEFFICIENTS FOR n,,^p 



Similarly, we will display the forms of Wilson coefficients emerged in the Eq. p9ti2"5|) after performing 
the Borel transformation. As been discussed in the text, only perturbative part contributes to the three- 
point function, 

= Fr'l + 0{as) + 0{l/mh). (Bl) 
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The relationship between Ff^^^, and p'?^^* are given as 



ppert 



(mc+mq)2 



(si -pI){s2-pI) 



dsi 



(B2) 



or 



r° ds2 ^«i7^e-^'/''^'li^e-"/*'^V'r*(«i,«2,g^) (B3) 



where the definition of Sj" is given in Eq. (|A4p . The manifest expressions of p'Y^^{i 
1, 4, 5, 6, 7, 9, 11, 12, 14, 15) are displayed as 
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{ml -w? + S2)(si + S2 - q^) - 2siS2 
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APPENDIX C: DECAY CONSTANTS OF J/i/' AND D^*] IN TWO-POINT QCD SUM RULES 



In this appendix, we would like to collect the sum rules for the decay constants of J and ^ for the 



comp 



as 



|45|, 



eteness of the paper. The decay constant of J /ip in the two-point QCD sum rules can be written 



58| 



i-2 2 - 



,0 

Aral 47r2 



4m 



CNl/2 



(1 + 



2mi 



)e m2 _|_ 



1 



1 



, a 



-G'^J0),(C1) 



s • • s ■ ■ 4M2 16 48 

where the nonrelativistic approximation for the gluon condensate has been adopted for the convenience of 
performing the Borel transformation. It is observed that the gluon condensate has tiny effect on the results 
of the form factors and hence are neglected in the sum rules of charmed mesons. The non-perturbative 
condensates used in the evaluation of the sum rues can be grouped as 



(C2) 



(0|^g|0) = -(1.65 ±0.15) X 10-2 GeV^ (g = u,(i), (0|ss|0) = (0.8 ± 0.1)(0|gg|0), 
(Ol^G^JO) = 0.005 ± 0.004GeV^ (0|g,a • Ggi|0) = mg(0|gigi|0), 

where = (0.8 it 0.2)GeV2 and the subscript "i" denotes the flavor of quarks. Based on the two- 
point sum rules of J/V' and the parameters showed above, we can derive the decay constant of J/ip 
as 337l:J^MeV, where we have combined the uncertainties from the variations of the Borel masses and 
threshold value for J/tp channel. 

The sum rules for the decay constant of Dq can be given by 
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(C3) 



)e I7^(0|^a-G(z|0), 



from which we can arrive at the decay constants of pseudoscalar charmed mesons as fo^ = 166l^QMeV 
and Jd, = 189l?oMeV. 

The decay constants of vector charmed mesons fo* in the framework of QCD sum rules can be 



calculated as 
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from which we can achieve the decay constants of vector charmed mesons as fo* = 240+|gMeV and 
fD* = 262t?2MeV. 
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